
Figures S1 A and S1B. UMAP of clusters highlighting macrophages (green, full line) and myofibroblasts (red, straited line). Cell-type clusters also include smooth muscle cells, fibroblasts, endothelial cells, granulocytes, NK/T cells, and B cells. (C) Breeding strategy used to create the macrophage-myofibroblast double-reporter (MAMY) mice. First, we crossed Postn MCM/MCM mice with Rosa26 tdTom a to/tdTom a to mice to obtain an established tamoxifen-inducible cardiac myofibroblasts lineage-tracer , . We then crossed these mice with a global monocyte/macrophage reporter-Cx3cr1 GFP/GFP . Triple-heterozygote mice were used as experimental animals (Postn MCM/+ Rosa26 TdTomto/+ Cx3cr1 GFP/+ ). (D) Experimental scheme whereby MAMY mice (Postn MCM/+ Rosa26 TdTomto/+ Cx3cr1 GFP/+ ) underwent sham or MI by permanent left anterior descending artery (LAD) ligation. To label myofibroblasts with tdTomato, mice were given daily tamoxifen injections following MI (red line). Cardiac function was tracked using 2D echocardiography at multiple time points following MI: 7 days prior to MI (baseline), and 7, 14, and 28 days post-MI (black arrows). Hearts were collected for histological fibrosis and immunofluorescence (IF) analysis on days: 2, 4, 7, 14, and 28 post-MI/sham (red arrows). (E) Ejection fraction (EF; %) analysis was performed on MI/Sham operated MAMY mice (Sham, n = 4; MI, n = 7). To extract the direct % change in EF between days 28 post-MI and baseline, we calculated the ΔEF value per animal . Red, MI; blue, sham. Results are represented as mean ± SEM. Statistical analysis performed using two-way ANOVA with Sidak’s adjusted p values (left) or two-tailed unpaired t test (right). Individual points represent individual biological replicates. (F) Representative histology images of MAMY hearts stained with sirius red on days 2, 4, 7, 14, and 28 post-MI/sham. Scale bar, 1 mm. Collagen, red; healthy myocardium, yellow. (G) MAMY mice fibrosis area/left ventricle (LV) in % was calculated per timepoint: Sham ( n = 12), day 2 ( n = 7), day 4 ( n = 4), day 7 ( n = 9), day 14 ( n = 4), and 28 ( n = 13) post-MI. Mix-max violin plot represent the median (middle full line) and quartiles (striated lines). Statistical analysis performed using one way ANOVA with Tukey’s multiple comparisons test. (H–K) Quantification of spatiotemporal distribution of cardiac-troponin-T cells (cardiomyocytes, cTnT+; purple), the lineage-traced (tdTomato+; red) myofibroblasts, and macrophages (GFP+, green) in Sham hearts and infarct zone of MI operated hearts. cTnT+, tdTomato+, and GFP+ cells were also measured at the suture area of 28 days post-MI hearts ( ,
Figure S1 C). (H) Representative immunofluorescence (IF) images of Sham ( n = 5; 1,487 ± 118.5 SEM cells per replicate) and MI operated hearts on days 2 ( n = 4; 1,088 ± 341.8 SEM cells per replicate), 4 ( n = 4; 1,027 ± 352.6 SEM cells per replicate), 7 ( n = 4; 1,792 ± 153.3 SEM cells per replicate), 14 ( n = 3; an average of 1,471 ± 275.2 SEM cells per replicate), and 28 ( n = 5; 1,627 ± 248.7 and 1,049 ± 315.9 SEM cells per replicate for infarct zone and suture, respectively) post-MI. Images are represented as either single channel for GFP, tdTomato, and cTnT, as a merged IF image and as a cell-type map that illustrates the identified cells as representative points. White scale bar, 100 μm. cTnT+ cells (I), GFP+ cells (J), and tdTomato+ cells (K) were quantified as count per 100 μm neighborhood or as % of cells per 100 μm neighborhood ( ,
Figure S1 C). Results are represented as mean ± SEM. In (I)–(K), infarct zone is denoted as a full circle, day 28 suture is denoted as an upside-down triangle. Statistical analysis performed using the Mann-Whitney test with Bonferroni adjusted p values (I–K). (L) Experimental scheme whereby adult Hsd:ICR (CD1) mice hearts underwent MI by permanent LAD ligation and further processed by flow cytometry . LV samples below the suture of MI or sham operated mice were collected at 6 different time points following injury at days: 0/Sham ( n = 11), 2 ( n = 5), 4 ( n = 5), 7 ( n = 7), 14 ( n = 5), and 65 ( n = 5). (M–Q) (M) Representative flow cytometry plots of the gating scheme used to identify (N) total cardiac macrophages (out of total immune cells- CD45 + ), (O) infiltrating cells (monocytes and neutrophils out of total immune cells-CD45 + ), (P) cardiac activated macrophages, measured by CD11b mean fluorescence intensity (MFI), and (Q) Resident macrophages (TIM4 + macrophages out of total macrophages. Data presented as mean ± SEM. Statistical analysis performed using one way ANOVA with Dunnett’s multiple comparisons procedure. " width="100%" height="100%">
Journal: Cell Systems
Article Title: Cold and hot fibrosis define clinically distinct cardiac pathologies
doi: 10.1016/j.cels.2025.101198
Figure Lengend Snippet: Acute myocardial infarction results in cold fibrosis (A) Mathematical model of the myofibroblast-macrophage cell circuit predicts three possible outcomes: hot fibrosis, cold fibrosis, and healing based on the abundance of these two populations (axes). The basin of attraction for the healing state is bounded by a separatrix. The scheme also denotes the growth factor signaling as either ON (black lines) or OFF (gray lines) at each specific fixed point. (B) Single-cell RNA sequencing (scRNA-seq) dynamic analysis of total macrophage and myofibroblast populations (% of total interstitial cells) in an adult mouse heart following-MI data obtained from Forte et al. Data are represented as cell composition fold change (FC) to day 0. These data are further described in Figures S1 A and S1B. UMAP of clusters highlighting macrophages (green, full line) and myofibroblasts (red, straited line). Cell-type clusters also include smooth muscle cells, fibroblasts, endothelial cells, granulocytes, NK/T cells, and B cells. (C) Breeding strategy used to create the macrophage-myofibroblast double-reporter (MAMY) mice. First, we crossed Postn MCM/MCM mice with Rosa26 tdTom a to/tdTom a to mice to obtain an established tamoxifen-inducible cardiac myofibroblasts lineage-tracer , . We then crossed these mice with a global monocyte/macrophage reporter-Cx3cr1 GFP/GFP . Triple-heterozygote mice were used as experimental animals (Postn MCM/+ Rosa26 TdTomto/+ Cx3cr1 GFP/+ ). (D) Experimental scheme whereby MAMY mice (Postn MCM/+ Rosa26 TdTomto/+ Cx3cr1 GFP/+ ) underwent sham or MI by permanent left anterior descending artery (LAD) ligation. To label myofibroblasts with tdTomato, mice were given daily tamoxifen injections following MI (red line). Cardiac function was tracked using 2D echocardiography at multiple time points following MI: 7 days prior to MI (baseline), and 7, 14, and 28 days post-MI (black arrows). Hearts were collected for histological fibrosis and immunofluorescence (IF) analysis on days: 2, 4, 7, 14, and 28 post-MI/sham (red arrows). (E) Ejection fraction (EF; %) analysis was performed on MI/Sham operated MAMY mice (Sham, n = 4; MI, n = 7). To extract the direct % change in EF between days 28 post-MI and baseline, we calculated the ΔEF value per animal . Red, MI; blue, sham. Results are represented as mean ± SEM. Statistical analysis performed using two-way ANOVA with Sidak’s adjusted p values (left) or two-tailed unpaired t test (right). Individual points represent individual biological replicates. (F) Representative histology images of MAMY hearts stained with sirius red on days 2, 4, 7, 14, and 28 post-MI/sham. Scale bar, 1 mm. Collagen, red; healthy myocardium, yellow. (G) MAMY mice fibrosis area/left ventricle (LV) in % was calculated per timepoint: Sham ( n = 12), day 2 ( n = 7), day 4 ( n = 4), day 7 ( n = 9), day 14 ( n = 4), and 28 ( n = 13) post-MI. Mix-max violin plot represent the median (middle full line) and quartiles (striated lines). Statistical analysis performed using one way ANOVA with Tukey’s multiple comparisons test. (H–K) Quantification of spatiotemporal distribution of cardiac-troponin-T cells (cardiomyocytes, cTnT+; purple), the lineage-traced (tdTomato+; red) myofibroblasts, and macrophages (GFP+, green) in Sham hearts and infarct zone of MI operated hearts. cTnT+, tdTomato+, and GFP+ cells were also measured at the suture area of 28 days post-MI hearts ( , Figure S1 C). (H) Representative immunofluorescence (IF) images of Sham ( n = 5; 1,487 ± 118.5 SEM cells per replicate) and MI operated hearts on days 2 ( n = 4; 1,088 ± 341.8 SEM cells per replicate), 4 ( n = 4; 1,027 ± 352.6 SEM cells per replicate), 7 ( n = 4; 1,792 ± 153.3 SEM cells per replicate), 14 ( n = 3; an average of 1,471 ± 275.2 SEM cells per replicate), and 28 ( n = 5; 1,627 ± 248.7 and 1,049 ± 315.9 SEM cells per replicate for infarct zone and suture, respectively) post-MI. Images are represented as either single channel for GFP, tdTomato, and cTnT, as a merged IF image and as a cell-type map that illustrates the identified cells as representative points. White scale bar, 100 μm. cTnT+ cells (I), GFP+ cells (J), and tdTomato+ cells (K) were quantified as count per 100 μm neighborhood or as % of cells per 100 μm neighborhood ( , Figure S1 C). Results are represented as mean ± SEM. In (I)–(K), infarct zone is denoted as a full circle, day 28 suture is denoted as an upside-down triangle. Statistical analysis performed using the Mann-Whitney test with Bonferroni adjusted p values (I–K). (L) Experimental scheme whereby adult Hsd:ICR (CD1) mice hearts underwent MI by permanent LAD ligation and further processed by flow cytometry . LV samples below the suture of MI or sham operated mice were collected at 6 different time points following injury at days: 0/Sham ( n = 11), 2 ( n = 5), 4 ( n = 5), 7 ( n = 7), 14 ( n = 5), and 65 ( n = 5). (M–Q) (M) Representative flow cytometry plots of the gating scheme used to identify (N) total cardiac macrophages (out of total immune cells- CD45 + ), (O) infiltrating cells (monocytes and neutrophils out of total immune cells-CD45 + ), (P) cardiac activated macrophages, measured by CD11b mean fluorescence intensity (MFI), and (Q) Resident macrophages (TIM4 + macrophages out of total macrophages. Data presented as mean ± SEM. Statistical analysis performed using one way ANOVA with Dunnett’s multiple comparisons procedure.
Article Snippet: For primary cardiac myofibroblast enriched cultures, cells were first isolated from adult (10-12 weeks) female ICR mice hearts using a neonatal dissociation kit (Miltenyi Biotec,130-098-373) and gentleMACS homogenizer (Miltenyi Biotec).
Techniques: RNA Sequencing, Ligation, Immunofluorescence, Two Tailed Test, Staining, MANN-WHITNEY, Flow Cytometry, Fluorescence

Figure S1 C). (J) Representative IF images of interstitial and perivascular sham ( n = 4; 1,833 ± 117.4 and 946.5 ± 71.68 SEM cells per replicate, respectively), and interstitial ( n = 6; 1,422 ± 156.7 SEM cells per replicate), replacement ( n = 6; 1,474 ± 286.5 SEM cells per replicate) and perivascular ( n = 6; 963.5 ± 68.77 SEM cells per replicate) fibrosis (fib) in hearts 28 days post-TAC. Images are represented as either single channel for: GFP (green), tdTomato (red), and cTnT (purple), as a merged IF image and as a cell-type map illustrated the identified cells as representative points. White scale bar, 100 μm. cTnT+ cells (K), GFP+ cells (L), and tdTomato+ cells (M) were quantified as count per 100 μm neighborhood or as % of cells per 100 μm neighborhood ( ,
Figure S1 C). Results are represented as mean ± SEM. Each dot represents a single biological replicate. Statistical analysis performed using the Mann-Whitney test with Bonferroni-adjusted p values when appropriate. " width="100%" height="100%">
Journal: Cell Systems
Article Title: Cold and hot fibrosis define clinically distinct cardiac pathologies
doi: 10.1016/j.cels.2025.101198
Figure Lengend Snippet: Chronic ventricular pressure overload results in hot fibrosis (A) Experimental scheme whereby MAMY mice (Postn MCM/+ Rosa26 TdTomto/+ Cx3cr1 GFP/+ ) underwent sham or transverse aortic constriction (TAC) to induce chronic pressure overload. To label myofibroblasts with tdTomato, mice were given 7 daily tamoxifen injections immediately following TAC, and on days 13–14 and 27–28 post-TAC (red line). Cardiac function was tracked using 2D echocardiography at baseline (−7) and on multiple time points following TAC: 7, 14, and 28 days (black arrows). Hearts were collected for histological fibrosis and immunofluorescence (IF) analysis 28 post-MI/sham (red arrows). Heart weight (HW) to body weight (BW) ratios were further analyzed 28 days following TAC (green arrow). (B–E) Temporal echocardiography measurements of MAMY mice following TAC/sham, including (B) aortic pressure gradient (mmHg), (C) left ventricle anterior wall diameter in diastole (LVAW;d, in mm), (D) left ventricle posterior wall diameter in diastole (LVPW;d, in mm), and (E) ejection fraction (EF; %). Results are represented as mean ± SEM. Statistical analysis performed using two-way ANOVA with Sidak’s adjusted p values. (F) HW/BW ratios (mg/gr) of MAMY mice sham ( n = 4) and TAC ( n = 9) hearts 28 days post-TAC. Mix-max violin plot represent the median (middle full line) and quartiles (striated lines). Statistical analysis performed using two-tailed unpaired t test. (G) Left: representative images of MAMY hearts sirius red stained sections, 28 days post-TAC/sham. Scale bar, 1 mm. Right: fibrosis in TAC hearts was quantified separately based on its histological presentation: interstitial, replacement and perivascular fibrosis. Scale bar, 100 μm. Collagen, red; healthy myocardium, yellow. (H) MAMY mice TAC interstitial ( n = 9) and replacement ( n = 8) fibrosis (fib) field of views (FOVs) were compared with healthy sham ( n = 4) myocardium (data presented as average %/FOV). Mix-max violin plot represent the median (middle full line) and quartiles (striated lines). Statistical analysis performed using one way ANOVA with Tukey’s multiple comparisons test. (I) MAMY mice TAC perivascular fibrosis ( n = 9) was compared with sham perivascular regions ( n = 4). Average fibrosis area/FOV (%) was calculated. Mix-max violin plot represent the median (middle full line) and quartiles (straited lines). Statistical analysis performed using two-tailed unpaired t test. (J–M) MAMY mice neighborhood quantification of cardiac-troponin-T cells (cardiomyocytes, cTnT+; purple), the lineage-traced (tdTomato+, red) myofibroblasts, and macrophages (GFP+, green) in TAC/Sham hearts ( , similar analysis demonstrated for MI hearts in Figure S1 C). (J) Representative IF images of interstitial and perivascular sham ( n = 4; 1,833 ± 117.4 and 946.5 ± 71.68 SEM cells per replicate, respectively), and interstitial ( n = 6; 1,422 ± 156.7 SEM cells per replicate), replacement ( n = 6; 1,474 ± 286.5 SEM cells per replicate) and perivascular ( n = 6; 963.5 ± 68.77 SEM cells per replicate) fibrosis (fib) in hearts 28 days post-TAC. Images are represented as either single channel for: GFP (green), tdTomato (red), and cTnT (purple), as a merged IF image and as a cell-type map illustrated the identified cells as representative points. White scale bar, 100 μm. cTnT+ cells (K), GFP+ cells (L), and tdTomato+ cells (M) were quantified as count per 100 μm neighborhood or as % of cells per 100 μm neighborhood ( , Figure S1 C). Results are represented as mean ± SEM. Each dot represents a single biological replicate. Statistical analysis performed using the Mann-Whitney test with Bonferroni-adjusted p values when appropriate.
Article Snippet: For primary cardiac myofibroblast enriched cultures, cells were first isolated from adult (10-12 weeks) female ICR mice hearts using a neonatal dissociation kit (Miltenyi Biotec,130-098-373) and gentleMACS homogenizer (Miltenyi Biotec).
Techniques: Immunofluorescence, Two Tailed Test, Staining, MANN-WHITNEY
![Cold fibrosis after MI is conserved in humans and in a clinically relevant porcine model (A and B) (A) Representative human left ventricle spatial transcriptomics slides (Visium) of patients following acute-myocardial infarction (MI) and non-transplanted donor hearts. Samples were divided based on time following MI as either: uninjured ( n = 10), early (days 0–15 post-MI; n = 6), and late (30 days+; n = 6). (B) Abundance of fibroblasts and myeloid cells was quantified based on deconvolution scores of cell types per spot . Myofibroblast were calculated as enrichment of the mean myofibroblast state score within spots with a minimal 10% value of cell-type abundance . Statistical analysis used Wilcoxon tests with Benjamini-Hochberg adjusted p values. (C) Experimental design of pig MI experiment: adult (3 months old) pigs underwent MI by temporarily occluding their LAD using a balloon catheter . Following reperfusion (balloon deflation), pigs were immediately treated with recombinant human Agrin (rhAgrin) or Saline control, in an antegrade trajectory. Injured pig hearts were collected at either day 3 ( n = 4 for rhAgrin; n = 3 for Saline) or day 28 ( n = 4 for rhAgrin; n = 3 for saline) and dissected to distinct tissue areas (infarct and remote zones). Remote and infarcted samples were subjected to bulk-mRNA sequencing and histology for fibrosis assessment. (D and E) Representative sirius red staining images are shown from (D) day 3 and (E) day 28 post-MI. Fibrosis was quantified as the average % fold change (FC) between infarct/remote zone sections for each pig individually. Fibrosis (day 3 or day 28) for Saline and rhAgrin hearts was measured using two-tailed unpaired t test. Scale bars: 1 mm. striated line denotes FC = 1. n.s, non-significant difference. (F and G) Heatmaps based on log 2 transformed normalized counts for all (upregulated and downregulated) differentially expressed genes (defined by |log 2 fold change| ≥ 1, p -adjusted value < 0.05 and max raw counts > 30) between remote and infarct zones for Saline and rhAgrin-treated hearts at day 3 (5,961 genes) (F) and 28 (1,079 genes) (G) post-MI. Rows represent genes and columns represent each biological sample and its spatial distribution according to infarct or remote zone. Data are represented as mean ± SD. (H) Hierarchical clustering per condition (day [3 or 28] +treatment [rhAgrin or Saline]), based on the 1,000 most variable genes. Triangles and circles represent remote and infarct zones, respectively. (I) Deconvolution of bulk-mRNA sequencing of pig hearts following MI of either rhAgrin (blue) or Saline-treated (red) samples. Macrophage and myofibroblast abundances were assessed by gene signatures as FC, between infarct and remote zones ( , C, and <xref ref-type=](https://pub-med-central-images-cdn.bioz.com/pub_med_central_ids_ending_with_2821/pmc11922821/pmc11922821__gr3.jpg)
Table S5 ). (J) Macrophage and myofibroblast abundances, based on deconvolution of bulk-mRNA sequencing (as in I). Treated samples were compared per time point (day 3 or 28), separately by two-tailed unpaired Student t test. Results are represented as mean ± SEM. " width="100%" height="100%">
Journal: Cell Systems
Article Title: Cold and hot fibrosis define clinically distinct cardiac pathologies
doi: 10.1016/j.cels.2025.101198
Figure Lengend Snippet: Cold fibrosis after MI is conserved in humans and in a clinically relevant porcine model (A and B) (A) Representative human left ventricle spatial transcriptomics slides (Visium) of patients following acute-myocardial infarction (MI) and non-transplanted donor hearts. Samples were divided based on time following MI as either: uninjured ( n = 10), early (days 0–15 post-MI; n = 6), and late (30 days+; n = 6). (B) Abundance of fibroblasts and myeloid cells was quantified based on deconvolution scores of cell types per spot . Myofibroblast were calculated as enrichment of the mean myofibroblast state score within spots with a minimal 10% value of cell-type abundance . Statistical analysis used Wilcoxon tests with Benjamini-Hochberg adjusted p values. (C) Experimental design of pig MI experiment: adult (3 months old) pigs underwent MI by temporarily occluding their LAD using a balloon catheter . Following reperfusion (balloon deflation), pigs were immediately treated with recombinant human Agrin (rhAgrin) or Saline control, in an antegrade trajectory. Injured pig hearts were collected at either day 3 ( n = 4 for rhAgrin; n = 3 for Saline) or day 28 ( n = 4 for rhAgrin; n = 3 for saline) and dissected to distinct tissue areas (infarct and remote zones). Remote and infarcted samples were subjected to bulk-mRNA sequencing and histology for fibrosis assessment. (D and E) Representative sirius red staining images are shown from (D) day 3 and (E) day 28 post-MI. Fibrosis was quantified as the average % fold change (FC) between infarct/remote zone sections for each pig individually. Fibrosis (day 3 or day 28) for Saline and rhAgrin hearts was measured using two-tailed unpaired t test. Scale bars: 1 mm. striated line denotes FC = 1. n.s, non-significant difference. (F and G) Heatmaps based on log 2 transformed normalized counts for all (upregulated and downregulated) differentially expressed genes (defined by |log 2 fold change| ≥ 1, p -adjusted value < 0.05 and max raw counts > 30) between remote and infarct zones for Saline and rhAgrin-treated hearts at day 3 (5,961 genes) (F) and 28 (1,079 genes) (G) post-MI. Rows represent genes and columns represent each biological sample and its spatial distribution according to infarct or remote zone. Data are represented as mean ± SD. (H) Hierarchical clustering per condition (day [3 or 28] +treatment [rhAgrin or Saline]), based on the 1,000 most variable genes. Triangles and circles represent remote and infarct zones, respectively. (I) Deconvolution of bulk-mRNA sequencing of pig hearts following MI of either rhAgrin (blue) or Saline-treated (red) samples. Macrophage and myofibroblast abundances were assessed by gene signatures as FC, between infarct and remote zones ( , C, and Table S5 ). (J) Macrophage and myofibroblast abundances, based on deconvolution of bulk-mRNA sequencing (as in I). Treated samples were compared per time point (day 3 or 28), separately by two-tailed unpaired Student t test. Results are represented as mean ± SEM.
Article Snippet: For primary cardiac myofibroblast enriched cultures, cells were first isolated from adult (10-12 weeks) female ICR mice hearts using a neonatal dissociation kit (Miltenyi Biotec,130-098-373) and gentleMACS homogenizer (Miltenyi Biotec).
Techniques: Recombinant, Saline, Control, Sequencing, Staining, Two Tailed Test, Transformation Assay